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Abstract 

The "anarchic" Randall-Sundrum model of flavor is a low energy solution to both the electroweak 
hierarchy and flavor problems. Such models have a warped, compact extra dimension with the 
standard model fermions and gauge bosons living in the bulk, and the Higgs living on or near 
the TeV brane. In this paper we consider bounds on these models set by lepton flavor violation 
constraints. We find that loop- induced decays of the form I — > I 7 are ultraviolet sensitive and 
uncalculable when the Higgs field is localized on a four-dimensional brane; this drawback does not 
occur when the Higgs field propagates in the full five-dimensional space-time. We find constraints 
at the few TeV level throughout the natural range of parameters, arising from u — e conversion 
in the presence of nuclei, rare \i decays, and rare r decays. A "tension" exists between loop- 
induced dipole decays such as u — ► and tree-level processes such as u — e conversion; they have 
opposite dependences on the five-dimensional Yukawa couplings, making it difficult to decouple 
flavor-violating effects. We emphasize the importance of the future experiments MEG and PRIME. 
These experiments will definitively test the Randall-Sundrum geometric origin of hierarchies in the 
lepton sector at the TeV-scale. 
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I. INTRODUCTION 



The Standard Model (SM) of particle physics is a remarkably successful description of 
nature. However, it contains several unsatisfactory features. In particular, there are many 
hierarchies built into the model that have no a priori explanation. The most famous of 
these is the huge separation between the electroweak and Planck scales. There have been 
many proposed solutions to this problem. One possibility is the Randall- Sundrum scenario 
(RS) ilj. In this model, our four- dimensional space-time is embedded into a five- dimensional 
anti de-Sitter space. The extra "warped" fifth dimension is compactified on an orbifold. This 
space-time is described by the metric 

ds 1 = e-^ri^dx^dx" - r 2 c d<j) 2 , (1) 

where — it < < vr. Three-branes are placed at the orbifold fixed points = and = it 
(and its reflection at = — 7r). The brane at = is called the Planck or ultraviolet (UV) 
brane, while the brane at = 7r is called the TeV or infrared (IR) brane. For sizes of the 
fifth dimension kr c ~ 11 — 12, the TeV scale is obtained from the fundamental Planck scale 
via an exponential warping induced by the anti de-Sitter geometry: M TeV = M p ie~ knrc . It 
was shown that this setup can be naturally stabilized 0]. The original model placed all SM 
fields on the IR brane. 

This scenario does not explain all unnatural parameters in the SM. The fermion Yukawa 
couplings, except for the top quark coupling, are small and hierarchical. The minimal RS 
model offers no solution to this flavor hierarchy problem. In addition, the flavor sector in 
the RS model is sensitive to ultraviolet physics, and requires a cut-off of roughly 10 3 TeV to 
avoid dangerous flavor- changing neutral currents. This is problematic, as the only cut-off 
available is the electroweak scale. 

One solution to this problem is to permit some or all of the SM fields to propagate in 
the full 5.D space jH 0, |^ . The only requirement for solving the gauge hierarchy problem 
is to have the Higgs field localized near the IR brane. This immediately presents a solution 
to the flavor hierarchy problem 0, H| , since the Yukawa couplings of the Higgs field to the 
fermions become dependent on the position of the fermion fields relative to the IR brane. By 
placing fermions at different positions in the 5D bulk, a hierarchy in the effective 4D Yukawa 
couplings can be generated even with anarchic 0(1) 5D couplings. These "anarchic" RS 
models set all diagonal and off-diagonal Yukawa couplings to 0(1). In addition, allowing 
fermions to propagate in the bulk suppresses the operators leading to dangerous flavor 
changing neutral currents 0|. Some collider |?J and flavor physics |H E| [HJ EH HI 
phenomenology of these models has been considered previously. 

Additional work is needed to make this scenario fully realistic. It was shown that the 
simplest formulation leads to large violations of the custodial symmetry in the SM [l3| . 
There are two known solutions to this problem. The first extends the bulk gauge symmetry 
to SU (2) l x SU (2) when broken by boundary conditions, a bulk custodial SU (2) symmetry 
is preserved [Til]. The second model introduces large brane kinetic terms to suppress precision 
electroweak constraints |l5|. Both solutions allow for the masses of the first Kaluza-Klein 
(KK) excitations to be as low as 3 TeV, generating interesting phenomenology which may 
be observable at the upcoming Large Hadron Collider (LHC). 

In this paper we probe the anarchic RS scenario by examining its effects on lepton flavor- 
violating observables. We study here a minimalistic model; we assume the SM gauge group, 
KK masses of a few TeV or larger, and an anarchic 5D Yukawa structure. We allow the 
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Higgs boson to propag ate in the full 5D space, which encompasses features found in several 
recent models [la, [l7[. Specific theories such as those mentioned above with a left-right 
symmetric bulk or large brane kinetic terms will predict slightly different effects than we 
find here, but we believe that our analysis captures the most important effects. We note that 
the flavor violation we study here is completely independent of neutrino physics parameters. 
We subject the anarchic RS picture to a complete set of experimental constraints: the rare /i 
decays \i — > and ^ — > e + e~e ± , the rare r decays r — > {e, fi} r y and tri-lepton decay modes, 
and fi — e conversion in the presence of nuclei. We find constraints on the KK scale of a few 
TeV throughout parameter space. Interestingly, there is a "tension" between dipole operator 
decays such as I — > I 7 and the remaining processes. They have different dependences on the 
5D Yukawa parameters, leading to strong constraints throughout parameter space. We also 
find that when the Higgs field is localized on the TeV brane, the dipole decays I — *■ I 7 are 
UV sensitive and uncalculable in the RS theory. This does not occur when the Higgs boson 
can propagate in the full 5 D sp ace-time. We emphasize the important role played by several 
future experiments: MEG [l8| , which will improve the constraints on [i — ► cy by two orders 
of magnitude; PRIME which will strengthen the bounds on fi — e conversion by several 
orders of magnitude; super-£> factories, which will improve the bounds on rare r decays by 
an order of magnitude. Measurements from these three experiments will definitively test 
the anarchic RS picture. 

We briefly compare our work to previous papers on lepton flavor violation in the RS 
framework. Reference [20] studied lepton flavor violation in a scenario where only a right- 
handed neutrino propagates in the full 5D spacetime. The studies in (§1 [T2| allowed all 
SM fermions and gauge bosons to propagate in the bulk. Reference (§[ did not incorporate 
custodial isospin, and therefore considered KK masses of 10 TeV, while the paper [12] 
considered a model with structure in the 5D masses and Yukawa couplings. None of these 
studies considered a bulk Higgs field. They also did not address the UV sensitivity of dipole 
decays in the brane Higgs field scenario, nor did they discuss the tension between tree-level 
and loop-induced processes. We also present a more detailed study of future experimental 
prospects than previous analyses. 

This paper is organized as follows. In Section 2 we present our notation and describe the 
model. We discuss in Section 3 the fi — e conversion and tri-lepton decay processes, which 
are mediated by tree-level gauge boson mixing. We discuss the loop-induced decays I — > I 7 
in Section 4. In Section 5 we present our Monte Carlo scan over the anarchic RS parameter 
space. We summarize and conclude in Section 6. 

II. NOTATION AND CONVENTIONS 

In this section we present our notation and describe the model we consider. The basic 
action is 



The Lagrangian for gauge fields in the bulk, £ gauge , has been studied in [3j. £f er mion was 
presented in 0, 0, 0| using an IR brane Higgs boson; we will review the relevant formulae 
and discuss the transition to a bulk Higgs below. Our setup of the bulk Higgs field will 
follow the discussion in ll 711 . 
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A. Brane Higgs field 



We begin by considering the case of the Higgs field localized on the IR brane. The 
Lagrangian in this case is 

£ffi gg s = [DpH^Hy - V(H) - ^Yukawa] fa ~ 7r) + 5(0 + 7r)] , (3) 

where is the covariant derivative. /^Yukawa describes the Yukawa interactions with the 
fermions. The Lagrangian for bulk fermions was derived in 0,0,0]; it takes the form 

Acrmion = ME%T A D M V ~ Sgn(0) k C ^ . (4) 

where is the inverse vielbein. This Lagrangian admits zero-mode solutions. The c$ 
parameters indicate where in the fifth dimension the zero-mode fermions are localized: either 
near the TeV brane (c < 1/2) or near the Planck brane (c > 1/2). The 4D Yukawa couplings 
of these fermions are exponentially sensitive to the c$ parameters. We perform the KK 
decomposition of the fermion field by splitting it into chiral components, \I/ = + ^ ' Ri 
yielding 

2kr c \4>\ 



The c dependence becomes part of the KK wavefunction f c); explicit formulas for 
these wavefunctions can be found in 0,0,0]. 

The SM contains two types of fermions, corresponding to singlets (S) and doublets (D) 
under SU(2)l- In the SM, we require that the S fermions are right-handed while the D 
fermions are left-handed. However, in five dimensions we must have both chiralities. To get 
a chiral zero-mode sector we use the orbifold parity of RS models. In particular, we choose 
(Sr, Dl) to be even under the orbifold parity (Neuman boundary conditions) and (Sl, Dr) 
to be odd (Dirichlet boundary conditions). The odd fields will not have zero modes, and 
the even zero modes will correspond to the SM fermions. We now group these fermions and 
their first KK modes into the vectors 

*i = {Df\Df\sf\ 

K = (Sf\sf\Di\ (6) 

where i is a flavor index (i = e,fi,r) and / = 1...9. We will show in a later section that 
higher KK modes have a negligible effect on our results. 
The fundamental 5D Yukawa interaction is 

^Yukawa = ~^~D l L HS R . (7) 



Using the vectors in Eq. |H1 and substituting in the KK expansion of Eq. yields 

^Yukawa = ^¥ L H^ R + h.C, (8) 

where 



A4D \adFr 
A = I F L X 4D F L X 4D F R I . (9) 
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Each internal block is a 3 x 3 matrix, with 



L,R 



( f (i) /AO) 

' Je L ,R/ Je L ,R 



V 



\ 

f W ,A0) o 

JHl,rIJvl,r u 

o / 



(10) 



These should be evaluated on the TeV brane, since that is where the Higgs is localized. We 
find 



A 



1 - 2q)(1 - 2cj) r i-(a+cj) v xij 



e f (0) f (0) . y 

W ' kr /i h 'I' - \j ( e l-2 Ci _ 1< _ ^ 



(11) 



where e = e 71 "^ and there is no sum over It is straightforward to write down the mass 
matrix for the fermions: 



M MqF r 
M = I F L M F L M F R M KK 
M KK 



(12) 



where Mq = ^A^ D is the zero mode mass matrix. Mkk is a diagonal matrix that contains 
the KK masses. Mq is not diagonal. We can diagonalize this zero mode mass matrix in the 
usual way, by constructing a biunitary transformation (Ul,Ur) so that Md = UlM U r is 
diagonal. We can embed this rotation into the full matrix above by multiplying on the left 
by diag^L, 1, 1) and on the right by diag(L r j ? , 1, 1). This gives 



M 




(13) 



f l u{m d . 



A factor of 



We have set^A R = U L M F R = M D U R F R and ^A L = F L M U R 
-7= was extracted to make it easier to match to the Yukawa matrix. Notice that the middle 
entry can also be written in terms of the diagonal zero-mode matrix: Ai = FlM F r = 
FjJJ\My)IJ R F R . From now on, we will use this expression. To find the Yukawa matrix A in 
this basis, we just divide Eq. [T^by and set M K k = 0. We note that this implies we are 
considering the exchange of a complex Higgs boson, which is equivalent to the exchange of 
the physical Higgs boson and the longitudinal component of the Z. The diagonalization of 
this mass matrix is discussed in the Appendix. 



B. Bulk Higgs field 

We now discuss the changes that occur when we allow the Higgs to propagate in the 
full 5.D space. A new coupling of the Higgs boson to the fermion KK states exists: 
HSl D R +h.c. This is not present in the brane Higgs case because the Sl and D R wave- 
functions vanish identically on the TeV brane due to the Dirichlet boundary conditions. The 
fermion mass matrix becomes 



M = A x M KK ■ (14) 
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Al,r,i are not the same as in the brane case; they now include overlap integrals of the 
KK and zero-mode fermion wavefunctions with the Higgs wavefunction. A2 represents the 
wavefunction overlaps between the first KK modes of the right-handed doublet and left- 
handed singlet leptons; the explicit expressions as well as the details of diagonalizing Ai can 
be found in the Appendix. We note that all of the A are proportional to the 4D Yukawa 
couplings. 

Our discussion of the bulk Higgs field will follow the presentation in (I?! . The 5D profile 
for the Higgs vev is 

Xh(<P) = N H e 2 °J v (ix T e kr ^) . (15) 

Here, xt is the solution of a root equation giving the tachyonic mass, Nh is a normalization 
factor, a = kr c (j), and v is the index of the solution. We will simplify this further for 
our discussion by using the asymptotic expansion of the Bessel function for large index, 
J u (z) ~ z v . Using this expansion gives the following normalized profile: 



This satisfies the constraint 



2 / dcjie^xM), (17) 







where the factor of 2 comes from the [— tt, 0] integration. In our analysis we will vary the 
index u, without worrying about its dependence on the model parameters in jl3|. This also 
makes a connection with the A 5 composite Higgs models in which is approximately 
realized in this framework as v = 0. We can also make a direct comparison to the TeV 
brane Higgs scenario, which is realized by v — > 00. 

We will now study the effect of the bulk Higgs field on the gauge boson sector. We begin 
with the action 



S gauge 



j d 5 xV^GG MN {D M H) ] D N H. (18) 



Performing a standard KK decomposition, and expanding H = vxh / V^c, we arrive at the 
mass matrix 



ml 



f a mn A^A^\ (19) 

with 



2 

m,n=0 



a mn = 4vr f T d<t>e- 2 °x 2 H X (m) X (n) - (20) 
Jo 

The x^ are the usual gauge wave-functions, which can be found in 0. We note that x^ — 
1 / \^2tt. We show in Fig. ^ the elements fi = a^i of this mixing matrix. The expectation is 
that as v — > 00, these should approach the brane Higgs values of (— l) l+1 \ / 27ikr c w ±8.42, 
assuming the value kr c = 11.27; this is indeed what occurs. 

We must now study the fermion sector, particularly what form the 4D Yukawa couplings 
take in terms of the 5D values. We begin with the action 



Sf, H = I d 5 * V=GrfMh (21) 
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FIG. 1: fi, the off-diagonal elements of the gauge boson mass matrix that describe the mixing of 
the zero-mode with the i-ih KK-mode. 

where the //-dependent prefactor is included to reproduce the correct 4D Yukawa coupling 
as v — ► oo. The zero-mode fermion wave-function is e 2<T / \Fr~ cl where 



/ 



(o) 



kr c (l - 2c) 

2 ( e {\-2c)kr c it _ \\ 



(22) 



Inserting this into the action, and expanding H as before, we find the following expression 
for the 4D Yukawa: 



A 



ID 



A 5 p (1 - 2c) 

e (l-2c)fer c 7r _ ^ 

A 5 p (1 - 2c) 

p (l-2c)fer c 7T _ 1 



y/kr c {l + v) 



'XHe 



-2ar 



l + V 



e (2+i/-2c)fcr c 7T _ i 



^/ e 2(l+v)kr c -K _ I 2 + V - 2C 



(23) 



For simplicity, we have only presented the diagonal Yukawa coupling. To reproduce the 
brane Higgs diagonal Yukawa coupling in Eq. [TTJ the bracketed integral should reduce to 



,(1— 2c)kr c Tr 



as v — > oo. It is simple to check that this occurs. 



C. The anarchic RS parameters 

We discuss here the parameters of the anarchic RS model and give their natural values. 
We first note that Eq. ^2 relates the diagonal 5D Yukawa couplings to the fermion c pa- 
rameters through the measured fermion masses. The off-diagonal entries are removed after 
diagonalization with Ul,r- The preferred size of the Yukawa couplings can be determined 
by demanding consistency with Z — > bb measurements and by the size of the top quark 
mass; this yields \ 5 o ~ 2 [14]. We assume three couplings (Y e ,Y^,Y T ) of this approximate 
magnitude. The size of these couplings implies c > 1/2 for all three leptons, indicating that 



they are localized near the Planck brane. For simplicity, we take cl = c R . We note that this 
range of c is the appropriate one for first and second generation quarks also; for the third 
generation, c& £ = q l ~ 0.45, while Cb R ~ 0.5 and c tR ~ |l4j . 

We can also estimate the natural sizes of the Ul jR matrix elements. For illustration, we 
consider here a two-family scenario; it is straightforward to extend this example to three 
families. Assuming an anarchic RS scenario, so that all of the X^ D ~ 0(1), we can use 
Eq. ^2 to write the 4D Yukawa matrix as 

where we have assumed for simplicity a symmetric 5D Yukawa matrix. Assuming 0(1) 
Yukawa couplings, the functional dependences of the fermion masses on the wave-functions 
are 

m e ~ /<°> 2 , m„ ~ 4 0)2 , (25) 
while the mixing matrices take the form 

U~\ I- (26) 



We therefore find that | C^x2 1 ~ \i- Including the r then gives \Uis\ ~ \P^- and \U> £ 
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The diagonal entries \Ua\ ~ 1. We will assume mixing matrix elements of these 
approximate magnitudes in our analysis. 



D. Operator Matching 

We discuss in this subsection the formalism we will use to compare the RS predictions to 
the experimental measurements. Our presentation closely follows the discussion in jij]. Tri- 
lepton decays of the form I — > £1^3 and fi — e conversion are mediated by tree-level mixing 
with heavy gauge bosons and generate four-fermion interactions, while I — > I 7 occurs via 
a loop-induced dipole operator. We can parameterize these effects in the following effective 
Lagrangian: 

-£ eff = C R (q 2 )^-e R a^F^ L + C L (q 2 )-^-e L a^F^ R 
AG F 

+7/|~ [93(e R ^fi R )(e R 'j^e R ) + 9i{e L ^ ^ L ){e L ^^e L ) 

+ 9s(eR r f/j, R )(e L 'Y ll e L ) + g 6 (e L 'j l "fx L )(e R 'y^e R )] + h.c. (27) 

The form factors 1 Ci yR (q 2 ) and the couplings ^ are then computed in a straightforward 
matching procedure. We will discuss this computation in detail in the following two sections. 



HSG Tfi 
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III. TRI-LEPTON DECAYS AND fi - e CONVERSION 



In this section and the next we study the predictions that the minimal RS model makes for 
lepton flavor violation. We focus on processes in the muon sector, such as fi~ —>■ e + e~e~ and 
li—e conversion in the presence of nuclei, and rare tau decays of the form r — > li^h currently 
being studied at BABAR and BELLE. The dipole-mediated decays will be discussed in the 
next section. 

The dominant effects arise from flavor non-diagonal couplings of the zero-mode Z-boson. 
Contributions from exchange of the Higgs boson are suppressed by small fermion masses, 
and we will show later that those coming from direct KK exchange are suppressed by a 
large fermion wave-function factor. There are also contributions to these processes from the 
dipole exchanges denoted by Cl,r in Eq. EH but these are loop-suppressed and small in the 
parameter space of interest. We also find that KK- fermion mixing effects are sub-dominant 
in the parameter space of interest. We derive here the relevant couplings. We denote the 
physical basis by Z ,Z 1 , and the gauge basis by Z^\Z^ l \ For simplicity, we restrict our 
discussion here to the first KK level; in our analysis we include the first several modes. After 
diagonalizing the gauge boson mass matrix, we find that these are related via 

Z^ = Z + f^ Zl , ZU = Zl -f^Z . (28) 

f parameterizes the mixing between the zero and first KK level. With a brane Higgs field, 
/ = y/2kirr c ~ O(10). A plot of / for a bulk Higgs field is shown in Fig. ^ The couplings 
between the zero-mode fermions and Z^ 1 ' are determined by the appropriate overlap integral. 
We define the ratio of these couplings to the SM couplings as a e , a M , and a T , where = 
ag SM ; the ol% are then given by 

a l = 2V2^ #e CT X (1) [if] 2 . (29) 
Jo 

Since the fermion wave-functions are localized at different points in the bulk, the a« differ, 
but they are all roughly O(0. 1) in magnitude. We present a plot of the ati in Fig. El In the 
fermion flavor basis, the matrix which describes the Z^ l > couplings takes the form 

/ a e \ / e F \ 
g SM (e F ,fi F ,f F ) a M Lu F . (30) 

\0 a T J \t f J 

We must first rotate the fermions to the mass basis. As was explained in the last section, we 
introduce unitary matrices Ul, Ur, so that L M = UlL f , R m = UrR f , where L F denotes the 
left-handed flavor basis-vector, Lm the left-handed mass basis-vector, etc. The flavor-basis 
coupling matrices C[ R = g^n diag(a e , a M , a T ) are rotated to Cl^r = Ul,r C[ r U* l r . The 
flavor-violating couplings are the off-diagonal entries of C^ ^; we find 

(l)/ie ( TT L,R TT L,R* . TT L,R TT L,R* . TT L,R TT L,R* \ 

9l,r = 9l,r U 2 { a e + U 2 j + C/ 13 ' a T J , 

9l'r = 9l,r \U 2 { U 3 { a e + U 22 ' U 32 ' + U 23 ' U 33 ' a T J , 

(l)re / tt L,RjtL,R* . TT L,R TT L,R* . TT L,R TT L,R* \ / 01 \ 

9l 'r = 9l,r \y x { U 3 { a e + U 12 U 32 + U 13 U 33 ' a T J , (31) 
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where g^R are the usual SM couplings. We can use the unitarity of Ul,r to rewrite these as 



„(!)*» 

9l,r 


— 9l,r 






- ot e ) 




[a T 


- 0t e ) 


9l,r 


— 9l,r 


ttL,RjtL,R* 

U 21 U 31 




~ «m) 


+ vtfvg* 




- <v) 


(l)re 

9l'r 


— 9l,r 






- oc e ) 




y 0l T 


- 0i e ) 



(32) 



Using Eq. |2HJ the couplings to Z are obtained via multiplication by —fm z /Mj CK : g^ R = 

—fm z /M] CK g ( i^ e , etc. The couplings to Z\ are identical to those in Eq. EH to leading order 
in the gauge boson mixing. 

We now use these to derive the flavor-violating couplings g 3 _ 6 of Eq. [2*71 



tie 
93 



tie 
94 



ue 
95 



fie 



2g R 

2g L 
2g R 



ae . . ^ 

9r +ot e g R 



{l)fie 



m 



2 1 



M 2 KK 



91 + a e9 L 



(l)fie 



rrir 



M 2 KK , 



9r + «eyR 



(33) 



These are for fi — e flavor violation; similar expressions hold for t — fi and r — e. The first 
term on each line is from the Z coupling, while the second is from direct Z\ exchange. 
Substituting in the expressions from Eq. "221 we find 



lie 
93 



"2<7 



rrir 



R 



M 2 KK 



(f - a e ) [U*U**( afl - a e ) + U*U%{a T - a t 



(34) 



and similar expressions for the other couplings. Since / ^> \a e \, we can neglect the direct 
KK exchange effect. 

We will study the decays fi~ — > e~e + e~ , r~ — > t~ — > e~e + e~ , r — ► /i~e + e~, and 

t — > e~fti + fU~. The remaining rare r decays studied at BABAR and BELLE, r — > e~fi + e~ 
and r — > fi~e + fi~ , require an additional flavor- violating coupling than those above, and are 
therefore highly suppressed. The relevant branching fractions from j2Jj are 



BR{p - 3e) = 2 (|C| 2 + \9T?) + ICI 2 + K 



/ie 1 2 



BRir^ Zp) = {2 (\g7\ 2 + \97\ 2 ) 



\9s 



rui2 



+ \g^\ 2 }BR(r^evv) 



BR(t - 3e) = {2 (|^ e | 2 + |yl e | 2 ) + |^ e | 2 + I^T} BR{t -> evv) 
BR(r^f,ee) = { |^T + \q7? + \g?\ 2 + ISeT} ^i?(r - evv), 
BRir^efifi) = {|<£ e | 2 + \g™\ 2 + + bST} BR{r - ew). 



(35) 



We have used the fact that BR(fi —>■ evv) = 1 in writing these expressions. The fi — e 
conversion rate is given by [22! ] 



2p e E e G 2 ,mfo 3 Z e 4 // Q 2 v 



(36) 
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where Gf is the Fermi constant, a is the QED coupling strength, and the remaining terms 
are atomic physics constants defined in 12 211 . Numerical values for titanium, for which the 
most sensitive limits have been obtained |23|. can be found in |2l| . 

We will present a detailed scan of the anarchic RS parameter space in a later section. For 
now, to provide some guidance as to what scales these rare decays can probe, we perform 
a few simple estimates. We set the 5-D fermion Yukawas to the values suggested by 5-D 
Yukawa anarchy, Y e = = Y T = 2. We also use the intuition described in the previous 
section to set the mixing matrix entries to the values 



U 



L,R 
11 



1, u 



L,R 
12 



m e ttL,R 
j ^13 



(37) 



and similarly for the remaining rows of Ul,r] for this estimate, we set the phases of these 
elements to zero. We choose a value of kr c = 11.27. We include the first 3 KK modes in 
this estimate, and we have checked that adding more does not affect our results. Employing 
these approximations, we check what limits can be obtained on Mkk from each process. We 



impose the following bounds: BR(fi — > 3e) < 10 -12 , which is the current PDG limit 24l |: 
B n nnv < 6.1 x 10~ 13 , which is the strongest constraint obtained by the experiment SINDRUM 
II |2i| . For the rare tau decays, we employ the strongest constraints from either B AB AR or 
BELLE, which are BR(r — > I1I2I3) < 2 x 10~ 7 for each mode [2^. We present the bounds 
on Mkk for both the brane Higgs model and the bulk Higgs scenario with v = in Table 
The limits from BR(fi — > 3e) and B conv already probe the multi-TeV region, similar to that 





Brane Higgs 


v = 


BR(n -> 3e) 


2.5 TeV 


2.0 TeV 


Bconv 


5.9 


4.7 


BR(t -» 3/i) 


0.40 


0.33 


BR(t — > iiee) 


0.36 


0.30 


BR(t -» 3e) 


0.10 


0.09 


BR(t — ► e/ifj,) 


0.09 


0.08 



TABLE I: Constraints on the first KK mode mass, Mkk, coming from various measurements for 
both a brane Higgs field and for the bulk Higgs case with v = 0. The bounds on Mkk are in TeV. 



possible at the LHC Although the limits from rare r-decays are lower, they probe different 
model parameters which describe the third generation. These bounds will also improve as 
the .B-factories acquire more data. We will show that these bounds are generic throughout 
the entire parameter space in a later section. 



IV. DIPOLE OPERATOR MEDIATED DECAYS 

We now compute the decays of the form I — > Z'7, which are induced at the loop level 
by the diagram shown in Fig. 121 For simplicity, we discuss the decay \i — > e^. It is simple 
to translate our expressions into results for r decays. The dominant contributions to these 
amplitudes come from exchange of a Higgs boson and KK fermions. This is because these 
diagrams contain terms proportional to the fourth power of the fermion wave-function ratio 
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H 



P 



V 



1 



FIG. 2: The Feynman diagram generating the dipole operator which mediates I — > Z 7 decays. Z* 
are the physical KK leptons. We have specialized to fi — > in the figure. There is a similar 
diagram with L «-> i?. 



f(l)/f(0) 



For c = 1/2, this ratio is f e il = 2irkr c w 70; it grows rapidly for c > 1/2, 



the values relevant for the muon and the electron. This strong dependence on the fermion 
wave-function was first noted in j2^| . There are also contributions coming from loops of KK 
Z bosons and KK fermions. However, as argued in reference for the case of the KK gluon 
contribution to radiative quark decays, the flavor structure of this diagram is approximately 
aligned with the 4D Yukawa matrix and hence gives a suppressed contribution. The KK 
fermion-Higgs diagrams have a different flavor structure than the 4D Yukawa matrix. 
The amplitude for the diagram in Fig. El is 



i 



d 4 k 



(2 



u(p)(iA e 0i 



V m KK 



V 2 - 



■{ik itJ p)u(p) 



k 2 — m 



2 

H 



-eA 



5^A e o, / 
4 J 



d A k 



(2vr) 4 (f. 



-A iM o 



u(p), (38) 



where = p^ + k and Ajj are the Yukawa matrices. We will assume the external lines 
are massless, which is valid up to subleading corrections in 1/ f e>fl . We have denoted the KK 
fermion masses by M%- K . At each KK level, there are two vector-like fermion pairs for each 

flavor with masses M$ K and M$ K , as is clear from Eq. El The splitting of these masses 
through mixing will be important in evaluating this contribution. It is straightforward to 
evaluate this integral to find 



A(jjl -> er, q 2 ) 



1 



2m, 



(39) 



where 
-*C(q 2 ) 



lem 



16tt 2 



2_j ^ e0i 



'0 jo 



X-z 



dz I dy 



M%{1 - z) 



q 2 y(l-y-z)-(l-z)M j 



kk zm 2 H 



Aj M o. 



(40) 
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We now set q 2 = to derive 



*^»)^E*.^v, (4i) 



KK 

where I(x) = 1 — x + (9(x 2 ). The branching fraction becomes [2 

12tt 2 

BQi - e T ) = — — [\C L (0)\ 2 + \C R (0)\ 2 ] , (42) 

where we have inserted the helicity labels L, R on C . These helicity labels dictate which 
elements of the Yukawa matrix A should be used; we will make this explicit in the following 
discussion. We now consider separately the brane and bulk Higgs field cases. We will find 
that the brane Higgs prediction for I — > I 7 is not calculable because it is sensitive to cut-off 
scale physics, while for the bulk Higgs case we can use our 5D effective field theory to make 
robust predictions. 



A. UV sensitivity for the case of brane Higgs field 

The leading contribution in Eq. El with m# = and I(x) = 1, vanishes up to factors 
suppressed by 1/f 2 for a brane Higgs field because of the Yukawa matrix structure. With 
rriH = 0, we are only considering contributions proportional to 1/M K ' K . This mass splitting 
is cancelled by shifts in the Yukawa couplings to all orders in v/M KK . The leading result 
therefore comes from 1/(M^'^) 3 contributions, and we must consider the m 2 H terms to 
obtain these. The diagonalization of the fermion mass matrix in Eq. El is discussed in the 
detail in the Appendix. The result of this analysis is the following mass splitting: 



1 3Ai 



This yields the following coefficients of the dipole operator: 



(43) 



The Yukawa structures entering Cx and Cr differ by a hermitian conjugate. 

However, it turns out that this result is masked by cut-off effects. A similar ultravio- 



let sensitivity of Higgs- fermion KK loops was also noted in The expected one- loop 

contribution from a given set of KK modes is finite with size 



C KK \ 2 1 



m* IQti 2 M 2 KK 



(45) 



For simplicity, we have not included the relevant mixing matrix elements in this estimate. 
Although the actual one-loop result for a brane Higgs field vanishes for m# = 0, we cannot 
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find a symmetry that requires this, and we expect it to be an accident of the one- loop result. 
The sum over two independent KK modes would have given a logarithmic divergence at 
one-loop: 

qkk 

' " oc log N KK ~ log (A 5D /k) ~ log [A 5D /M KK j . (46) 



ml 



Here, Nkk is the total number of KK modes in the 5D effective theory, A 5 £> is the 5D cut-off 
of order 10 19 GeV, and A 5 £> is the warped-down 5D cut-off of order TeV. Similarly, M^k 
is roughly the warped-down curvature scale k. To obtain this logarithmic divergence, it is 
crucial that KK fermion-Higgs couplings in the sum are independent of the KK index. We 
expect that higher-loop contributions are strongly power divergent because of the increasing 
number of sums over KK modes, and are as important as the one-loop result provided the 
cut-off scale physics is strongly coupled. 

This divergence structure can be more easily seen using power-counting in the 5D theory. 
Since the 5D Yukawa coupling has mass dimension [X^n/k] = [— 1], the loop expansion for 
fi — > e7 has the form 



r KK 1 / \ 



m 2 16tt 2 \M KK 



log 1 k ) 167T 2 p A5D + 



(47) 



In this expression, we have replaced the scale k ~ 10 18 GeV by its warped-down value M^k 
in the overall coefficient. By simple dimensional analysis, the one-loop contribution can be 
log-divergent and the two-loop contribution is quadratically divergent; in KK language, the 
power divergence at two loops can be seen from the independent sums over 4 KK modes. 
The two-loop result is comparable to the one-loop prediction if the cut-off physics is strongly- 
coupled: A\ D jk 2 x \\ D I (167r 2 ) ~ 1. Therefore, the KK loop contribution is not calculable 
in this case. 

Based on the above discussion, we also expect the higher- dimensional operators in the 
5D theory coming from physics at the cut-off scale to be important. The relation between 
the warped-down 5D cut-off in the Yukawa sector and the KK scale for a brane Higgs field 
is A 5D ~ MkkJ (4vt/A 5 £)) , based on power counting of the 5D loop factor. To obtain the 
cut-off operator, we replace M KK in Eq. ESI by the cut-off scale A 5D , and the loop factor 
by ~ 1, since the cut-off effect has no loop suppression. This shows that /i — > e'j is an 
UV sensitive observable for a Higgs field on the TeV brane. We can only parameterize the 
contribution as: 

2 

CZr = x U^f, (48) 

where a is an unknown, 0(1) coefficient, and we have included the appropriate mixing matrix 
element. 

We now show that we can reliably calculate dipole induced decays for a bulk Higgs field. 
The Yukawa coupling in this case has mass dimension [X^o/Vk] = [— 1/2], so the loop 
expansion is instead 



^L,R J- / *5D 



ml 16vr 2 \VM^ 



1 , 1 Agp /A^\ (J_\ Xj D - 
M KK 167r 2 M KK \ k ) + Mft K 



From this 5D power-counting, we see that the one-loop KK contribution is finite. The two- 
loop result is logarithmically divergent, but is smaller than the one-loop prediction by ~ 0.1 
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provided X^d ~ 4. Three-loop and higher contributions are power-divergent and comparable 
to the two- loop result, but are again smaller than the one-loop effect. 

Thus, in the bulk Higgs case, the KK effect is calculable. The effects from cut-off scale 
operators are suppressed, and we can reliably make a prediction using the RS theory. In our 
numerical analysis, we will include dipole decays for the bulk Higgs field case. For the brane 
Higgs scenario we will simply neglect them, since we cannot make a reliable prediction. 



B. Contributions from a bulk Higgs field 



We now consider the scenario when the Higgs boson is allowed to propagate in the bulk. In 
this case, the KK mode result is not overwhelmed by cut-off scale operators. The = 
limit does not vanish for a bulk Higgs. We make this approximation in our discussion, 
since the corrections are 0{m 2 H / M KK ). We first work out the Yukawa structure appearing 
in Eq. |^ Using the results in the Appendix for the two KK fermions appearing in the 
diagram of Fig. |2l we find 
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(49) 



i/j, 



In this expression we must sum over I = e, /i, r. To simplify this we use the splitting between 
the KK fermion masses derived in the Appendix: 
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KK 1V± KK 

We find the following results for the dipole operator coefficients: 

3em 



(50) 



Cl(0) 
C r (0) 



32% 2 M KK 



3em, 



32% 2 M KK 



[A iJ A 2 A i ; 



[A fi A 2 A L 



t 

eft 



(51) 



We note that in the limit of the Higgs boson being localized on the TeV brane, A 2 — > 0; the 
result vanishes in this limit, as required. 

An identical analysis can be performed for r — > /ry and r — > ej. We simply replace 
— > m T , change the indices of the Yukawa structure appropriately in Eq.^Hl and normalize 
the expression to the decay r — > evv. We now perform an estimate of the bounds similar to 
that performed in the brane Higgs case. We set Y e = Y iL = Y T = 2, and set the mixing matrix 
elements to their canonical values as described before. We also set v = 0. We impose the 
following bounds on each of the three dipole decays: BR(n — > ej) < 1.2 x 10~ n , as obtained 
from |24j | ; B R(r — > fij) < 9 x 10 -8 , the stronger of the bounds coming from BABAR and 
BELLE [27]; BR{r — > ej) < 1.1 x 10~ 7 , again the stronger of the bounds coming from 
BABAR and BELLE |28|. We find the following constraints for the canonical parameters: 



(52) 



BR(n - 


-> ej) 


M KK 


> 15.8 TeV 


BR(t - 


-> ery) 


M KK 


> 0.9 TeV; 


BR(t - 


■> fij) 


M KK 


> 1.6 TeV. 
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The constraints, particularly from BR(fi — > cy), are quite strong. This arises in part from 
the large value of the Yukawa coupling, Y = 2, as we now discuss. 

C. Tension between tree-level and loop-induced processes 

We now discuss a tension between processes caused by tree-level gauge boson mixing 
such as n — e conversion and I — > hhh, and dipole operator decays. These have opposite 
dependences on the 5D Yukawa couplings, leading to strong constraints for all parameter 
choices. We first give a very simple scaling argument to motivate this, and then present 
numerical proof. 

Our scaling argument uses the dependence of each process on the zero-mode fermion 
wave-function evaluated at the TeV brane. We will work for simplicity in the large v 
limit, which mimics a brane-localized Higgs field. From Eqs. |221andEHl we find that the 
wave-function scales roughly as /j ~ 1/ x/Ajsd- The wave-function has weak c-dependent 
factors which we will ignore in this argument. The quantity that governs the flavor violation 
in gauge boson mixing is the difference between o^'s, as is clear from Eq.|^] In the definition 
of ai in Eq. |2H we can divide the overlap integral into two regions, one near the Planck brane 
and the other near the TeV brane, to show that the former is c-independent and that the 

latter carries the c-dependence and must be ai\ non - universa i ~ [fi WreVbrane ~ 1/%d- We 
therefore expect the non-universal part of ai, and hence the flavor violation, to decrease for 
larger Yukawa couplings, which is indeed what we observe in Fig. EJ For the dipole mediated 
decays, recall that in Section |IV] we claimed that the operator coefficients Cl,r scaled as 
Cl,r ~ V[//°' ) ] 4 ~ ^5D) this can be verified using Eq. loTl and the results in the Appendix. 
The constraints coming from I — > I 7 decays will increase with larger Yukawa couplings, the 
opposite dependence of the tree-level processes. 

0. 1 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 



0.0 — 




FIG. 3: The ratios of the zero-mode fermion couplings to over their SM values, for x = e, pb, r, 
as functions of the Yukawa couplings Y x . 

To exhibit this behavior we present in Table HP the bounds on the first KK mode mass for 
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canonical mixing angles, v — 0, and for the two choices of Yukawa strength Y e = Y il = Y T = 
1,2. We show the two most constraining processes, fi — e conversion and BR(fi — > e^y). The 
dependence on the Yukawa couplings agrees with our simple estimate above. We will find in 
the next section that this leads to strong constraint throughout the entire model parameter 
space. 





Y = 1 


Y = 2 




6.7 TeV 


4.7 TeV 


BR(fi -► erf) 


8.0 


15.8 



TABLE II: Constraints on the first KK mode mass, Mkk> coming from \x — e conversion and 
BR(fi — > e*y), for canonical mixing angles, v = 0, and for Y = 1, 2. 



V. MONTE-CARLO SCAN OF THE ANARCHIC RS PARAMETER SPACE 

In this section we present our Monte-Carlo scan of the RS parameter space, to determine 
in detail how well the RS geometric origin of flavor can be tested by current and future 
lepton flavor-violation experiments. 

We first describe the ranges over which we scan the various RS parameters. The scenario 
introduced in the previous sections contains the following free parameters: Y e , Y^, Y T , the 
overall Yukawa couplings for the electron, muon, and tau; U^' , the elements of both the left 
and right-handed mixing matrices; the KK mass M^k- We make the following assumptions 
in our scan. 

• We restrict the Yukawa couplings to the range Y x e [|,4]. As discussed before, the 
natural value is Y x ps 2. Values larger than 4 begin to invalidate the perturbative 
expansion, while values smaller than 1/2 introduce an unnatural hierarchy in the 
model. We explained in the previous section that flavor violation cannot be removed 
by making the Yukawa couplings either large or small, due to tension between tree-level 
and loop-induced processes. 

L R 

• We implement the anarchy of 5-D couplings in our scan, which indicates that U u ' ~ 1, 
U 12 ' ~ ^m e /m^, C/ 13 ' ~ ^/m e /m T , etc. We fix U u ' =1, and define the canonical 
values 

ut, = M, u; 3 = M, u S , = J^. (53) 

y m M y m T Y m r 

We then vary = (3f 2 U^ 2 , with /3f 2 G [1/4,4]. We independently vary U^ R , and 
U^ R in a similar fashion. Again, we restrict the values to these ranges to insure no 
unnatural hierarchies in model parameters. We generate phases for the six independent 
jjl,r j n ran g e 2tt]. 

L R 

• We approximately implement unitarity of the mixing matrices by setting U 2 { = 

~ (Ul2 R )\ etc - This 

assures that unitarity is maintained up to corrections of the 
level ^Jrne/m^, ^Jm^/m T , which is sufficient for the scan performed here. 
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We scan over the following fifteen independent parameters: the three Y x , and the six complex 
mixing matrix elements U^, U^ R , and We generate 1000 sets of fifteen random 

numbers, and distribute them in the ranges indicated above for fixed Mkk- We perform 
two separate scans, one for a brane Higgs field and one for a bulk Higgs with v — 0. The v 
dependence of the bulk Higgs field bounds is studied separately. 



A. Scan for the brane Higgs field scenario 

We first perform a Monte-Carlo scan of the parameter space of the brane Higgs scenario. 
As discussed in Section I1V1 dipole decays of the form I — I 7 are UV sensitive. We do not 
consider these decays in the brane Higgs case, which leaves us with /1 — e conversion, /i — > 3e, 
and t — > li^h- 

We first study the muonic processes /a — > 3e and fi — e conversion. We show in Fig. 0] 
scatter plots of the predictions for BR(fi — > 3e) and B conv coming from our scan of the 
RS parameter space, for the A" A' scales Mkk = 3, 5, 10 TeV. The most sensitive probe is 
the SINDRUM II limit of B conv < 6.1 x 10~ 13 0]. This rules out a large fraction of the 
parameter space for Mkk < 5 TeV, and restricts the allowed parameters even at 10 TeV. 
The PDG limit of BR(fi — > 3e) < 10~ 12 is less severe: although it rules out a large fraction 
of the M K k = 3 TeV parameter space, most of the M K k = 5 TeV space is still allowed. We 
note there is an almost perfect correlation between the RS predictions for the two processes. 
This is not surprising; it is clear from Eqs. ESI and EHl that they depend almost identically 
on the same mixing angles. 

This result has implications for both the aesthetic appeal of the anarchic RS flavor picture, 
and the observation of this physics at the LHC Although points with Mkk < 3 TeV are 
still allowed, it is clear from Fig. 0] that the model as formulated in our scan prefers KK 
masses of 5 TeV or larger. Increasing the KK scale to these higher values introduces a large 
fine-tuning in the electroweak symmetry breaking sector and is therefore not favored 
With such large KK masses, many associated states will also be too heavy to observe at the 
LHC. The other method of avoiding these constraints, reducing the Uh' R matrix elements 
to the appropriate level, implies either some additional structure or fine-tuning in the 5-D 
Yukawa matrix. We have studied the minimal model here, and it seems likely that more 
structure in the 5-D Yukawa matrix is needed for a completely natural description of the 
first and second generation flavor pattern in the brane Higgs case. 

Another sector of the RS flavor picture to explore is that involving the third generation r. 
This tests different model parameters than the muonic processes. We show in Fig.EJa scatter 
plot of the RS predictions for BR{r — > 3e) and BR{r — > 3/i) for M K k = 1 TeV, together 
with the best limits coming from BABAR and BELLE. The lowest AA-scale allowed by 
electroweak precision tests in anarchic RS models is typically a few TeV. The 5-factories 
are beginning to probe this region in the mode r — > 3/i. There are plans to build a super- .B 
factory with an integrated luminosity approaching 10 ab _1 2^|. The projected limits from 



this experiment are included in Fig. |SJ Both the r — > 3/i and r — > 3e modes at a super- .B 
factory will constrain the anarchic RS parameter space. The LHC also has sensitivity to rare 
r decays |3(|; however, the projected sensitivities are slightly weaker than the current B- 
factory constraints, and have not been included. The expected sensitivities to rare r decays 
at a future linear collider are also weaker than the limits set by the S-factories. Although 
the Mkk ~ 1 TeV scales probed with r — > /1M3 decays are lower than those constrained by 
/i — e conversion and \i 3e, we stress that different model parameters are tested by each 
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FIG. 4: Scan of the \i — > 3e and \i — e conversion predictions for Mkk = 3, 5, 10 TeV. The solid 
and dashed lines are the PDG and SINDRUM II limits, respectively. 



set of processes. 



B. Scan for the bulk Higgs field scenario 



We now present the results of our scan over the bulk Higgs parameter space. For the scan 
we set v = 0, which mimics the composite (or A 5 ) Higgs model of [l6(; we present separately 
the v dependence of the most important constraints. 

We again begin by considering muon initiated processes. The constraints from /j —>■ 3e 
and /i — e conversion are highly correlated, as we saw in the previous subsection. Since 
the bounds from \x — e conversion are stronger, we focus on this and fi ej. We show in 
Fig. scatter plots of the predictions for BR(fi — > ej) and B conv coming from our scan of 
the RS parameter space, for the KK scales Mkk = 3, 5, 10 TeV. For /x — > we include 
both the current constraint from the Particle Data Group 24]| and the projected sensitivity 
of MEG [IH- The current bounds from \i — > ej are quite strong; from the Mkk — 3 TeV 
plot in Fig. El we see that only one parameter choice satisfies the BR{ji — > ej) bound. 
This point does not satisfy the [l — e conversion constraint. We can estimate that it would 
satisfy both bounds for Mkk > 3.1 TeV. In our scan over 1000 sets of model parameters 
the absolute lowest scale allowed is thus slightly larger than 3 TeV. Also, a large portion of 
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Mifif = 1 TeV 




FIG. 5: Scan of the r — ► 3e and r — > 3/i predictions for Mkk = 1 TeV. The solid and dashed lines 
are the current .B-factory and projected super-i? factory limits, respectively. 

the parameter set at both 5 and 10 TeV conflict with these bounds. We again find the need 
for a KK scale of Mkk > 5 TeV or additional structure in the mixing between the first and 
second generations to satisfy the experimental constraints for a significant fraction of model 
parameter space. In Fig. 0we present the anarchic RS predictions for r — > /ry and r — > ej, 
together with current and future 5-factory constraints, for Mkk = 3 TeV. The r — > /ry 
mode is currently probing the few TeV range, while r — > will begin to test the anarchic 
RS scenario during the running of a super- £? factory. 

To study the sensitivity of the bulk Higgs field scenario to the location of the Higgs boson 
in the fifth dimension, we show in Fig. |H] the dependence of the \l — e conversion rate and 
BR(fi — > ej) on v. We set the mixing angles to their canonical values, and show results for 
Y x — 1, 2 and Mkk = 5, 10 TeV. The fi — e conversion results are weakest for v = 0, and 
quickly asymptote to the brane Higgs result as v becomes large. The variation of fi — > 
with v is more intricate. The vanishing of the calculable component of this process as the 
Higgs boson is moved towards the TeV brane, discussed in Section IIVI is clearly seen in 
Fig. |H1 However, we expect cut-off effects to become more important for large v. There 
is a strong dependence of the process on the position of the Higgs field for small u, with 
the result varying by an order of magnitude for < v < 5. The v — case is again the 
most favorable choice. Since UV sensitivity of the model is reduced for a bulk Higgs field, 
and since the experimental constraints are weakest for v = 0, we conclude that there is a 
preference for models of the type presented in [lif . 

C. Future sensitivities of MEG and PRIME 

Finally, we emphasize here the importance of future searches for /i — e conversion by 
PRIME and \i — > &-/ by MEG. Our analysis has shown that with some small tuning of 
parameters, particularly for those describing the mixing of the first and second generation, 
KK scales of 3 TeV are allowed by current measurements. Alternatively, KK scales of 5 
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FIG. 6: Scan of the /i — ► and /i— e conversion predictions for Mkk = 3, 5, 10 TeV and v = 0. The 
solid line denotes the PDG bound on BR(fi — > ej), while the dashed lines indicate the SINDRUM 
II limit on fj, — e conversion and the projected MEG sensitivity to BR{[i — > ey). 



TeV are permitted with completely natural parameters. Super- 5 factory searches for rare 
r decays will not significantly constrain scales Mkk > 5 TeV. The LHC search reach for 
the new states predicted by the anarchic RS scenario is expected to be around 5-6 TeV. It 
is therefore difficult to definitively test the RS geometric origin of flavor using data from 
-B-factories and the LHC. 

Searches for \i — e conversion and \x — > are already starting to require slight tunings of 
the model parameters. The limit on BR(fi — > ej) is projected to improve from 1.2 x 10~ n 
to 10~ 13 after MEG, while the constraint on /i — e conversion is projected to improve to 
10~ 18 after PRIME. The bounds on Mkk that these constraints lead to are shown in Fig. El 
We have plotted the projected bounds as a function of the overall scale of the mixing angles; 
we have set = n^nie/m^, U^ R = n^m e /m T , etc., and have varied k in the range 

[0.01,1]. This tests how far from the natural parameters these experiments will probe. We 
observe that MEG will probe Mkk — 5 TeV down to mixing angles 1/10 times their natural 
sizes. PRIME will test Mkk — 20 TeV down to mixing angles 1/10 times their natural sizes, 
and will probe Mkk < 10 TeV down to mixing angles 1/100 times their canonical values. 
Together, these experiments will definitively test the anarchic RS explanation of the flavor 
sector. 
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FIG. 7: Scan of the r — > /ry and r — ► e7 predictions for Mkk = 3 TeV and i/ = 0. The solid and 
dashed lines are the current S-factory and projected super- B factory limits, respectively. 
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FIG. 8: v dependence of the RS predictions for \i — e conversion and /x — > for canonical mixing 
angles and for several choices of Y x and Mkk- In the right panel, the Y = 1, Mkk = 5 TeV and 
the y = 2, M Xi f = 10 TeV lines overlap. 



VI. SUMMARY AND CONCLUSIONS 

In this paper, we have studied lepton flavor violation with the SM propagating in a 
warped extra dimension. The principal motivation for this model is a solution to the Planck- 
electroweak hierarchy problem. Interestingly, there is also a solution to the flavor hierarchy 
of the SM. The large differences in the quark and lepton masses and mixing angles can 
be explained by differing profiles of SM fermions in the extra dimension, even though the 
5D Yukawa coupling are of the same size without any structure. These profiles can vary 
substantially with small changes in the 5D fermion masses; no large hierarchies are required 
to account for the flavor hierarchy in the SM. Since the Higgs field is localized near the TeV 
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FIG. 9: Projected bounds on Mkk coming from MEG (left) and PRIME (right) for v = 0. We 
have set the mixing angles to k times their canonical values, and have varied k in the range [0.01, 1] 
for Y x = 1,2. 



brane, the small masses of the first and second generations are explained by their localization 
near the Planck brane. 

The localization of fermion fields at different points in the extra dimension leads to flavor 
violation upon rotation to the fermion mass basis. The assumption of anarchic 5D Yukawa 
couplings implies that the mixing angles are related to the ratios of fermion masses. We can 
therefore estimate the leptonic mixing angles without a model of neutrino masses, unlike 
in the SM. The flavor violating couplings are proportional to the 4D Yukawa interactions. 
Therefore there is an analog of the GIM mechanism in the anarchic RS picture. However, 
the sensitivities of lepton flavor violating experiments are large, so we expect significant 
constraints. Bounds from electroweak precision measurements currently constrain the KK 
scale to be Mkk > 3 TeV, approximately. 

To derive the implications of lepton flavor violating measurements for the anarchic RS 
scenario, we perform a Monte Carlo scan over the natural parameter space of this model: 
0(1) Yukawa couplings and 0(1) variations of the mixing angles around their predicted size. 
We study both the case where the Higgs boson is localized in the TeV brane and when it is 
allowed to propagate in the full 5D spacetime. We study the processes \i — > 3e, r — > hhh, 
fi — e conversion, and dipole decays of the form I — > I 7. In the brane Higgs case, cut-off 
effects render the dipole decays uncalculable in the 5D RS theory; this arises from the fact 
that the 5D Yukawa couplings in this case have mass dimension [—1], and cut-off scale 
effects are as large as those from KK modes. The bulk Higgs case does not suffer from this 
drawback. 

We find strong constraints throughout the entire natural RS parameter space. The mini- 
mal allowed KK scale is 3 TeV, and this is permitted only for a very few points in our scan. 
In the bulk Higgs case, this occurs partially because of a tension between the tree-level me- 
diated fi — e conversion process and the loop-induced decay \i — > ej. These processes have 
opposite dependences on the 5D Yukawa couplings, making it difficult to decouple the effects 
of flavor violation. There are a couple of possible ways to avoid these constraints. First, 
the KK scale can be raised slightly to 5 TeV, which allows large regions of the natural RS 
parameter space to be realized. However, this increases the fine-tuning in the electroweak 
sector, and will make it difficult to find the KK states present in this model at the LHC. 
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Another possibility is to reduce the leptonic mixing angles slightly, implying some struc- 
ture in the 5D Yukawa matrix and indicating that the observed flavor structure cannot be 
generated completely via geometry. 

There are also several possible model-building possibilities to relax these constraints. 
Models with custodial isospin based on the gauge structure SU(2) L x SU(2) R x U{V)b-l 
contain an additional Z and possibly additional fermions. The coupling of the Z to the 
SM fermions is model-dependent j3l]], and can possibly be used to cancel some of the flavor- 
violating contributions we have studied. These models also contain an additional right- 
handed neutrino that contributes to loop-induced dipole decays. There is no zero-mode 
partner of this right-handed neutrino, and this contribution is therefore independent of the 
neutrino mixing parameters. Even an 0(1) suppression suffices to reduce the KK scale to 
the 3 TeV level, opening up more parameter space for study at the LHC. 

The definitive test of whether the observed flavor structure can be explained by the 
anarchic RS scenario will come from future lepton flavor violating measurements. 5-factories 
are currently probing mixing in the third generation using rare r decays. These constraints 
will improve by an order of magnitude with data from a super- B factory, probing KK scales 
up to 5 TeV. These measurements probe different model parameters than fi—e conversion and 
rare /x decays, and are therefore complimentary to these other experiments. Improvements 
in the sensitivities of /x — > and /x — e conversion of several orders of magnitude will 
be accomplished by the future experiments MEG and PRIME, respectively. They will 
definitively test the geometric origin of flavor structure; for example, PRIME will probe 
KK scales of Mkk > 10 TeV down to model parameters 1/100 of their natural size. These 
experiments will either confirm or completely invalidate this geometric origin of flavor. 

In conclusion, the anarchic RS picture is an attractive solution to both the electroweak 
and flavor hierarchies in the SM. Measurements at the LHC, at future 5-factories, and with 
the experiments MEG and PRIME will determine whether it is indeed realized in nature. 
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APPENDIX 

In this Appendix we present the expressions that appear in the fermion mass and Yukawa 
coupling matrices. We focus on the case of a bulk Higgs field; the brane Higgs mass matrix 
can be obtained by taking the appropriate limits, as discussed below. 

The fermion mass matrix for the bulk Higgs scenario is given by 



M 



( M ° ^ 
V M KK 



\ 

M KK 
A 2 



(54) 
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Mb is a 3 x 3 diagonal matrix containing the masses of the zero-mode leptons, and Mkk is 
the diagonal matrix with the KK masses before mixing. We assume that it is proportional 
to the identity matrix, as the deviation from this limit is small. The other entries can be 
expressed using the following two overlap integrals: 



(55) 
(56) 



) ^3 1 



Xh is the Higgs vev profile given in Eq. [THl while the f n ((f>; cj) can be found in 0]. The 
remaining matrices in Eq. |^] are 



V2 

v 

71 



A% = U Llk (U{M D U R ) k > x F k \ 



^2 



(U{M D U R y k U ] Rkj x F k \ 



A? = (U{M D U R f x Gl 
{U ] L M D U R f x G% 



(57) 
(58) 

(59) 
(60) 



where there is no sum over the indices i,j but there is over the index k. 

To see how this reduces to the brane Higgs case, we replace the Higgs wavefunction 
with a delta function on the TeV brane. This sets G R = via the boundary conditions 
/!(</> = tt) = 0, so that A 2 = as in Eq. ED Also, F lJ = F 3 \ since the i flavor cancels out 
of the ratio in Eq. and G l [ = F l F- \ again matching our results for the brane Higgs. 

We now discuss the diagonalization of this matrix and the fermion Yukawa coupling 
matrix. We first diagonalize the lower 2x2 block containing the KK masses, to remove 
the mixing between the KK fermions. We then diagonalize the full 3x3 matrix, to remove 
mixing between zero and KK modes. We include only the leading v/Mkk corrections. 

We first consider the following simple 2x2 matrix, which simulates the lower block of 
Eq.El 

T = (Vy)- < 61) 

We assume x, y < 1. This matrix is diagonalized by the following unitary transformation: 



V 



y/2 VI 



4 



"(1 + ^ 



(62) 



VTV^ is diagonal with eigenvalues ±(1± ^(x + y)) to leading order in x and y. We now make 
this the lower 2x2 block of a diagonalization matrix V, identifying X = x — y - Al_A - 
and compute Md = VMV^\ 



Mkk 



M 



D 



( 



M, 



D 



n/2 R 



V2 



1 + 



M KK + ^. 




4 - 







x - 

4 - 



-M KK + 



(63) 
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We now diagonalize the zero-KK mixing. We accomplish this with the following unitary 
transformation matrix: 



Y 



I 1 

v/V2 Al 



v/V2 
M K k 



A 



75' 



M KK 



R 72 
1 





;i-r) 



v/V2 
M KK 



A 



R 72 


1 



/ 



(64) 



where T = Al+ ^^ 4M ° . This removes the off-diagonal elements of the mass matrix. The 
eigenvalues are shifted at 0(v/Mkk)', this will be important for the Yukawa matrix below. 
The full diagonalization matrix is YV. We must also include the phase rotation P = 
diag(l, 1,-1) to make the eigenvalues positive. We can then determine the masses of the 
KK fermions to first order in v/Mkk'- 



M 



M 



(i) 

KK 

(2) 
KK 



M KK + 



M KK - 



Ai + A 2 



Ai + A 2 



(65) 
(66) 



These expressions for the KK masses are used when computing the amplitude for fx — > ery. 
The expression is valid for both the brane and bulk Higgs scenarios. 

We now determine the Yukawa coupling matrix. The Yukawa matrix A in the flavor basis 
is obtained by dividing by -^g and setting M KK = in Eq. Multiplying YV K{YV)^ P , 
we obtain the Yukawa matrix in the mass basis 



/ 



A 



D 



A 



ID 



1+4 



X _ A 2 -A/ 
4 M KK 



1-4 



A 2 -A/q 
M KK 



1 

\ 73 



1+1$ 



-Mq 



Mkk 



| _ [£_ A 2 -M 
1 4 M KK 



A, 
A, 



We do not include the lower 2x2 block since we do not need it here. Note that A 2 
so we can drop the dependence on the zero mode mass matrix in the off-diagonal 
These correspond to subleading contributions suppressed by l// 2 . 



(67) 
>M , 
terms. 
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